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Links in a practical network may have different functions, which makes the original network a combination 
of some functional subnetworks. Here, by a model of coupled oscillators, we investigate how such functional 
subnetworks are evolved and developed according to the network structure and dynamics. In particular, we 
study the case of evolutionary clustered networks in which the function of each link (either attractive or repulsive 
coupling) is updated by the local dynamics. It is found that, during the process of system evolution, the network 
is gradually stabilized into a particular form in which the attractive (repulsive) subnetwork consists only the 
intralinks (interlinks). Based on the properties of subnetwork evolution, we also propose a new algorithm for 
network partition which is distinguished by the convenient operation and fast computing speed. 

PACS numbers: 89.75.Hc, 05.45.Xt 



The past decade has witnessed the blooming of network sci- 
ence, in which one important issue is to explore the interplay 
between the network structure and dynamics 111 01 • While the 
influences of the network structure on dynamics have been in- 
tensively studied in the past [31, recently attentions have also 
been paid to the influences of the network dynamics on struc- 
tureji.e. the evolution of complex networks driven by dynam- 
ics ||l 1 H 1 11 El . In Ref. H it has been shown that, 
rewiring network links according to the node synchroniza- 
tion, a random network can be gradually developed to a small- 
world network. In Ref. JtII it has been shown that, driven by 
node synchronization, the weight of the network links can be 
developed to a particular form in favor of global network syn- 
chronization. Besides network evolution, dynamics has been 
also used for network detection, e.g., detecting the modular 
structures in clustered complex networks S 0. HI IH [Till ■ 

It has been well recognized that links in a practical network 
are usually different from each other In previous studies, this 
has been mainly reflected in the variation of the weight of the 
network links, i.e., the weighted network fill. Weighted net- 
work, however, describes only the case of single-function net- 
works, i.e. all links in the network have the same function, but 
failing to describe the situation of multi-function networks in 
which the network links have the diverse functions. A type 
of commonly seen multi-function networks in practice is the 
cooperation-competition network (CCN) lll2i [isll . in which 
the network links are divided into two groups of opposite 
functions. For instance, in the nervous network of the human 
brain, the synapses are roughly divided into two groups, exci- 
tatory and inhibitory, which play the contrary roles to the neu- 
ron activities 1IT2I1 . Another typical example of CCN is the re- 
lationship network shown in the prisoner's dilemma game, in 
which each suspect may either cooperate with (remain silent) 
or defect from (betray) the other suspects lUsll . 
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For multi-function networks like CCN, to facilitate the 
analysis, it will be convenient if we treat the different groups 
of links separately. That is, we pick out links serving the same 
function and, together with their associated nodes, construct a 
small single-function network. In this way, a multi-function 
network can be decomposed into a number of functional sub- 
networks, while each supports a unique function to the system 
behaviors. Here an interesting question is: How do these func- 
tional subnetworks co-evolve with each other and develope 
into their "adult" forms according to the system properties, 
e.g., the network structure and dynamics? 

To mimic the evolution of the functional subnetworks, we 
propose the following model of coupled phase oscillators, 

N 

Oi tJi+e^ay[sin(6'j-6'i)ey+cos(6'j-6'i)(l-ey)]. (1) 

Here, i,j = l,2,...,Ai" are the node indices, e is the uniform 
coupling strength. 9i and LOi are the instant phase and intrin- 
sic frequency of the ith oscillator, respectively. The network 
structure is represented by the adjacency matrix A = {fly }, 
in which = 1 if nodes i and j are directly connected, and 
aij — otherwise. E{t) = {ey(t)} is a time-dependent bi- 
nary matrix whose elements are defined as follows. Let 
be the instant phase of the local order parameter defined by 
the equation Ill4jl 

N 

r,e*'^'(*) =^a„e'«^W. (2) 

We set eij{t) = 1 if the difference between ijji and ijjj is 
smaller than a threshold D, otherwise we set [t] = 0. 

Different from the traditional models of coupled phase os- 
cillators, in Eq. ([T]l the coupling term is made up of two 
parts of the opposite functions. While the attractive coupling. 
Ha ~ sm{dj — 9i), is going to synchronize the connected 
nodes, the repulsive coupling, Hji — cos{9j ~ 9i), will work 
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against this tendency. These opposite functions, however, can- 
not coexist. That is, at any time instant each Hnk can only take 
on one type of coupHng function, either attractive (e^ — 1) or 
repulsive (e^ — 0). The attractive links, together with their 
associated nodes, constitute the attractive subnetwork, which 
is represented by the matrix B = A o ( " o " is the entry- wise 
product and / is the identity matrix). Similarly, we can con- 
struct the repulsive subnetwork, and represent it by the matrix 
R — A o {I — E). Because eijit) is being updated with the 
system dynamics, the two subnetworks, therefore, will also 
be changing with time. It should be noted that, despite the 
evolution of the subnetworks, the global network structure is 
kept unchanged, i.e., B{t) + R{t) = A. Imagine a complex 
network that is weakly coupled and there is no synchroniza- 
tion between any pair of nodes. It can be expected that, as 
the system evolves, the two subnetworks will be continuously 
updated in a random fashion. The question we are interested 
here is: What happens to the evolution of the subnetworks if 
the coupling strength is stronger? 

We start our investigation by considering the evolution of 
clustered networks (CN) [1]. A typical model of CN is the 
ad hoc network introduced in Ref. 1 15], which consists of 4 
clusters, each contains 32 nodes. In this model, each node on 
average has (fc) — 16 links, among which (fc;) links are con- 
nected to nodes within the same cluster, i.e. the intralinks, and 
(fcp) links are connected to nodes from different clusters, i.e, 
the interlinks. Since we are interested in the case of strongly 
coupled clustered networks, we use, without loss of general- 
ity, in our simulations the parameters e = 5 and (fcp) — 1. 
Meanwhile, to generate the matrix E{t), we use the threshold 
D — 0.45. (The influences of these parameters to the evo- 
lution will be discussed later) The natural frequencies and 
initial conditions of the oscillators are chosen randomly from 
the ranges [0, 1] and [0, 2tt), respectively. To monitor the evo- 
lution, we keep a record of the instant states of the oscilla- 
tors {{Oi{ty\) and the instant subnetwork matrices {B{t) and 

The evolution of the structures of the subnetworks can be 
described as follows. At the beginning, the two subnetworks 
have similar configuration, i.e. both are an abbreviated ver- 
sion of the original network [Figs. 1(a) and (d)]. This is 
because in a very short time the oscillators have not reached 
any synchronization, and therefore the matrices B and R are 
mainly determined by the initial conditions of the oscillators. 
But, due to the small value of D, the repulsive subnetwork has 
more links than the attractive subnetwork. Then, as time in- 
creases, the interlinks are gradually excluded from the attrac- 
tive subnetwork; meanwhile, the intralinks are excluded from 
the repulsive subnetwork [Figs. 1(b) and (e)]. The separa- 
tion of the two subnetworks, however, is not an even process, 
as some links may jump between the subnetworks repeatedly 
before settling down. Finally, at the time about T = 50, the 
subnetworks are stabilized into fixed structures and the evo- 
lution is complete. In this final stationary state, all intralinks 
(intralinks) of the network are included in the attractive (re- 
pulsive) subnetwork [Figs. 1(c) and (f)]. 

With the evolution of the subnetworks, the system dynam- 
ics is also changed. As shown in Fig. 2(a), with the increase 
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FIG. 1 : The evolution of functional subnetworks in the ad hoc net- 
work. (a)-(c) The evolution of the attractive matrix B{t). (d)-(f) The 
evolution of the repulsive matrix R{t). From time t « 50, all in- 
tralinks (interlinks) are contained in the attractive (repulsive) subnet- 
work, and the structures of the subnetworks will be stabilized. Nodes 
are rearranged according to the topological clusters. 
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FIG. 2: (Color online) For the same network as in Fig. 1, (a) the time 
evolution of the oscillator states and (b) a snapshot of the oscillator 
states at time t = 100. 



of time, the oscillators are gradually organized into 4 syn- 
chronous clusters. The pattern of the synchronous clusters 
is more evident in Fig. 2(b), where a snapshot of the oscil- 
lator states is taken at time t ~ 100. After a check of the 
node indices of the pattern, it is found that the organization of 
the dynamical clusters obey precisely the topological clusters. 
Moreover, in forming the synchronous clusters, it is found that 
the overlapping nodes, i,e. nodes which have interlinks, are 
more difficult to be synchronized than those internal nodes. 
This is clearly evidenced in the synchronizing process of the 
first cluster, where the few overlapping nodes are traveling 
among the synchronous clusters for an extra period before set- 
tling down (the black curves in Fig. 1(a)). Another interesting 
finding is that, in the stabilized pattern, the states of the syn- 
chronous clusters are well separated from each other, which, 
of course, is attributed to the repulsive coupling on the inter- 
links. 

A prominent feature of our evolutionary model is that the 
link functions are updated only by the local network informa- 
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tion, i.e., the phase of the local order parameter. While it has 
been commonly believed that the identification of the link at- 
tribute, i.e. intralink or interlink, relies on the knowledge of 
the global network information, e.g., the betweenness central- 
ity of the network links, it is somewhat surprising to see that 
here the identification can be accomplished by only the local 
network information. This interesting phenomenon can be ex- 
plained by a local mean-field theory, as follows. Let be an 
intralink of a clustered network. Since nodes inside a clus- 
ter are densely connected, nodes i and j thus are surrounded 
by a similar set of neighboring nodes. Because of the large 
overlap of their neighboring sets, the average phases ipi and 
il>j will have small difference, leading to the attractive cou- 
pling on the intralink, i.e., = 1. In contrast, if is an 
interlink, the nodes i and j will be surrounded by very differ- 
ent neighboring sets, which will generate a larger difference 
in the average phase, finally leading to the repulsive coupling 
on the interlink. 

We next discuss the influences of the network structure on 
the evolution. Having understood the critical role of the over- 
lapping neighbors in the evolution, we are able to predict that 
the above phenomena of subnetwork formation can be ob- 
served in any network of clear modular structures. To ver- 
ify this, we have studied the evolution of the other two typical 
network models. The first one is the overlapping clustered net- 
work studied in Ref. jlOtl . in which two larger clusters (each 
has 96 nodes) are mediated by two smaller complete clusters 
(each has 4 nodes). The smaller clusters have no direct con- 
nection, but each is connected to the two larger clusters by an 
equal number of links. By analyzing their neighbor sets, the 
network nodes are immediately classified into four groups: 
two for the larger clusters and two for the smaller clusters. 
Correspondingly, the oscillators are expected to be synchro- 
nized into 4 dynamical clusters. This is indeed what we have 
observed in the simulations [Fig. 3(a) and (b)]. The second 
model we have simulated is an ER network 1 1]. Since an ER 
network has no topological cluster, the neighboring sets of the 
network nodes thus are different from each other. Accord- 
ing to the neighbor-set analysis, this will lead to the repulsive 
couplings on the links, and generating the turbulent system 
dynamics. This is indeed what we have found at the beginning 
of the evolution [Fig. 3(c) and (d)]. The repulsive network and 
turbulent dynamics, however, are unstable. As shown in Fig. 
3, after a transient period, the repulsive couplings are quickly 
switched to the attractive couplings and the turbulent state is 
changed to the state of global synchronization [Figs. 3(c) and 
(d)]. The switching of the link functions and system dynamics 
suggest the dual properties of the ER network, i.e., it can be 
regarded either as containing no module or as containing one 
unified module. 

We go on to study the influences of other system param- 
eters on the evolution, including the coupling strength e, the 
threshold D, and the local dynamics. The numerical results 
show that, given the network has a clear modular structure, 
the link attribute can always be successfully detected by the 
local dynamics, despite the changes of e and D. Specifically, 
for the ad hoc network of Fig. 1, the system will always de- 
velop to the same functional subnetworks [Fig. 1] and dy- 
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FIG. 3: (Color online) For the overlapped clustered network, (a) the 
evolution of the system dynamics and (b) the stabilized link func- 
tions. In (b), the intralinks of the larger and smaller clusters are 
marked by black squares and red dots, respectively. The interlinks 
are marked by blue triangles. For a ER network of size A'^ = 128 and 
average degree (k) = 16, (c) the evolution of the system dynamics 
and (d) the stabilized link functions. The switching of the link func- 
tions and system dynamics are started at about t = 30. In (d), all 
network links are granted with the attractive coupling function. 



namical pattern [Fig. 2] in the parameter space constructed 
by £ G [0.5,8] and £) G [0.45, 7r/2]. Furthermore, the main 
feature of the evolution is independent of the specific form of 
the local dynamics, as has been verified by other nonlinear 
oscillators, such as the Logistic map and Lorenz oscillator. 
However, it should be pointed out that, by changing these pa- 
rameters, the transient process of the network evolution could 
be strongly affected, e.g., the transient states of the evolution. 

Finally, we discuss the possible application of the evolu- 
tionary model in network partition III Hill, [S [Hi ■ In 
network partition, the performance of an algorithm is mainly 
measured by the following three factors: ease of implemen- 
tation, accurate detection and fast computation. All these 
factors are well met in the evolutionary subnetwork model 
(ESM). Firstly, on the aspect of ease of implementation, the 
model employs only the local network information, and the 
partition is accomplished automatically by the system dynam- 
ics. In particular, at the end of the evolution, based on the 
states of the synchronous clusters, the topological clusters 
can be readily identified. Automatic detection and local net- 
work information are the main features of the ESM algorithm, 
which are also the major difference to the other dynamics- 
based algorithms |8, 9, 17]. Secondly, the ESM algorithm is 
fast in computing speed. The computational cost of the ESM 
method is estimated to be 0{NT), with N the network size 
and T the transient time of the evolution. To estimate the com- 
putational cost further, it is necessary to characterize the rela- 
tionship between T and N. Numerically, we have checked 
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FIG. 4: (Color online) For the ad hoc network of Fig. 1, (a) the 
transient time of the whole network, T, and the synchronizing time 
of the largest cluster in the nework, T', as a function of the network 
size, and (b) the variation of the mutual information (see Ref. 1 20] 
for details about mutual information) as a function of {kp) generated 
by the partition algorithms "Potts" 1 17], "ESM", "GKM" |9], "OCR" 
[8J, and "CNM" [18J. The coappearance matrix generated by ESM 
for (c) the Zachary's karate club and (d) the football club. The over- 
lapping nodes are marked by the squares. In (c) and (d), each data is 
averaged over 100 realizations. 



this relationship by increasing the size of the ad hoc network 
in the following two approaches. In the first approach, the 
number of the clusters are kept unchanged, but the size of the 
4th cluster is gradually increased from 32 to 864. In the sec- 
ond approach, the size of each cluster is kept unchanged, but 
the number of the clusters is increased from 4 to 30. The 
variations of T as a function of N are plotted in Fig. 4(a), 



together with the synchronizing time, T' , of the largest clus- 
ter in the network. Very interestingly, it is found that for both 
approaches we have T oc T'. That is, the transient time of 
the whole network is proportional to that of the largest cluster. 
Particularly, for the first approach we even have T ~ T' . Pre- 
vious studies have indicated that T' is mainly determined by 
the intrinsic properties of the cluster, e.g. the cluster size, in- 
stead of the global network properties ll2lll . This implies that 
T' ^ Mjnax, with Mjnax the size of the largest cluster in the 
network. Therefore, the computational cost of ESM is esti- 
mated to be proportional to 0{N M„iax)- Since for practical 
networks we generally have M,nax ^ N, the computational 
cost of ESM thus is estimated to increase linearly with the 
network size. Finally, on the aspect of detecting accuracy, the 
ESM algorithm works very well for clustered networks and 
reasonably well for fuzzy networks [Fig. 4(b)]. 

As applications of the new algorithm, we have tested the 
partitions of two empirical clustered networks. The first one 
is the Zachary's karate club network 1*22], which contains 34 
nodes and 78 links. The numerical result is plotted in Fig. 
4(c), in which the network is clearly divided into 4 clusters. 
Moreover, the overlapping node, i.e. the 10th node in the 
network, is also well characterized. These results coincide 
with that of Ref. |22|] . The second example we have tested is 
the football network fisl], which contains 115 nodes (teams) 
and 613 (matches) links. The numerical result is plotted in 
Fig. 4(d), in which the network is clearly divided into a num- 
ber of conferences (clusters). Again, the independent teams, 
which have equal number of games (links) with multiple con- 
ferences, are well identified by the overlapping nodes. These 
results coincide with that of Ref. ifTsI] . 

In summary, we have studied the evolution of functional 
subnetworks in clustered networks and proposed a new algo- 
rithm for partitioning networks. Hopefully, the finding that the 
network function is jointly determined by the network struc- 
ture and system dynamics could be helpful to the study of 
complex behaviors in functional networks. 
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